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Abstract 

Motivated by recent developments on random polymer models we propose a 
generalisation of reflected Brownian motion (RBM) in a polyhedral domain. This 
process is obtained by replacing the singular drift on the boundary by a continuous 
one which depends, via a potential U, on the position of the process relative to the 
domain. It was shown by Harrison and Williams (1987) that RBM in a polyhe- 
dral domain has an invariant measure in product form if a certain skew-symmetry 
condition holds. We show that (modulo technical assumptions) the generalised 
RBM has an invariant measure in product form if (and essentially only if) the same 
skew-symmetry condition holds, independent of the choice of potential. The invari- 
ant measure of course does depend on the potential. Examples include TASEP-like 
particle systems, generalisations of Brownian motion with rank-dependent drift and 
diffusions connected to the generalised Pitman transform. 

1 Introduction 

There has been much recent development on the study of random polymer models in 
1 + 1 dimensions [3,5,22,24,31,32]. These turn out to be closely related to random 
matrix theory on the one hand and the study of the KPZ equation [18], which was 
proposed to describe a class of surface growth models, on the other. An important 
role is played by an exactly solvable semi-discrete directed polymer model which was 
introduced in [25] and further studied in [3,21,22,32,33]. This polymer model can also 
be viewed as a network of generalised Brownian queues in tandem (or, equivalently, a 
TASEP-like interacting particle system). A important role is played by an analogue of 
the output or Burke theorem, which leads to a product-form invariant measure for the 
series of queues. 

Queueing networks [6, 13, 14, 28] provide examples of reflected Brownian motion 
(RBM) in a polyhedral domain, introduced and studied by Harrison and Williams 
[16,36]. The analogue of the output theorem in this more general setting is the exis- 
tence of an invariant measure in product form. The main result in [36] is that RBM in 
a polyhedral domain has an invariant measure in product form if and only if a certain 
skew- symmetry condition holds. 

Motivated by these facts, we introduce a multidimensional diffusion we call gener- 
alised RBM ( GRBM). Rather than giving the Brownian motion a singular drift whenever 
it hits one of the boundaries, we now impose a continuous drift. Its magnitude depends, 
via a potential U, on the position of the process relative to the polyhedral domain. A 
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special case, the exponentially RBM, is given by the choice U(x) = —e~ x , which corre- 
sponds to the generalised Brownian queue. We show that for the GRBM existence of 
an invariant measure in product form is still equivalent to the skew-symmetry condition 
of Harrison- Williams, irrespective of the function U. 

By introducing a parameter j3 (which can be viewed as inverse temperature), taking 
U(x) = -/3 _1 e _/3:E and letting (3 — > oo, we recover the RBM studied by Harrison- 
Williams. In this sense the GRBM process really is a generalisation of reflected Brow- 
nian motion. 

An example of GRBM is given by an analogue of the totally asymmetric simple 
exclusion process (TASEP) with a fixed number of particles. In the exponential case 
this corresponds to the above-mentioned generalised Brownian queues in tandem. The 
existence of an invariant measure in product form allows us to compute the asymptotic 
speed from equilibrium for this particle system, as well as the speed of the system 
when started from '— oo' in reverse order (the analogue of the so-called 'step initial 
condition'). In the exponential case this particle system can be viewed as a projection 
of a higher dimensional diffusion on generalised Gelfand-Tsetlin patterns which was 
introduced in [22] and shown to be closely related to the quantum Toda lattice. This 
higher dimensional diffusion in fact also fits into the framework of GRBM [23]. 

We also consider generalisations of Brownian motion with rank- dependent drift, stud- 
ied by Pal-Pitman [26] and draw connections to the 5-Bose gas recently studied and 
related to reflected Brownian motion in a Weyl chamber by Prolhac-Spohn [27]. 

This paper is organised as follows. In Section 2 we review reflected Brownian motion 
in a polyhedral domain. Section 3 is devoted to an analytic result that corresponds to 
the basic adjoint relations. The GRBM process is introduced in Section 4, where we 
also state the main probabilistic results. In Section 5 we discuss particular cases and 
examples. Section 6 contains proofs. 

2 Reflected Processes in a Polyhedral Domain 

We recall here the definition and main properties of the reflected Brownian motion 
in a polyhedral domain and define its exponential analogue. The former were first 
introduced and studied by Harrison and Williams [16,36]. A related process was 
studied even earlier by Harrison-Reiman [15], where the domain is always given by 
an orthant, but the driving Brownian motion is allowed to have a general covariance. 
Our generalisations cover both settings. 

2.1 RBM in a Polyhedral Domain 

Let us first discuss RBM in a general domain. In the Harrison- Williams setting 
the polyhedral domain G C M. d in which the process runs is the intersection of k > d 
half-spaces. More precisely let ri\,...,n^ £ R d be unit vectors and 6 £ 1* then the 
domain G is defined by 



k 



k 




We assume that G is non-empty and that each of the faces 



Fj = {x £ G : nj ■ x = bj 



} 
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has dimension d— 1. In general G may be bounded or unbounded, but we assume that 
{n\, . . . , rik} spans M. d , which means that no line can lie entirely within G. 

The reflections are defined by vectors qi, . . . , 6 ~R d such that qj ■ rij = for all j. 
We denote by N and Q the k x d matrices whose j th rows are rij and respectively. The 
requirement that {m, . . . , n^} spans all of M d is equivalent to existence of an invertible 
d x d submatrix N of N. 




Figure 2.1: Reflection vectors in a wedge 



Informally, reflected Brownian motion w in G may be described as follows: inside the 
domain G the process oj behaves like a standard Brownian motion with drift — [i £ M. d , 
at the boundary it receives an instantaneous singular drift pointing towards the interior 
- in direction Vj := qj + rij whenever it hits the face Fj - and it almost surely never 
hits any point in the intersection of two or more faces. 

In general such a process may not exist. The boundary of the state space is not 
smooth, and the directions of reflection are discontinuous across non-smooth parts of the 
boundary, so this does not fit within the Stroock-Varadhan theory [34] of multidimen- 
sional diffusions. Williams [36] showed that if the input data satisfy the skew- symmetry 
condition 

rij ■ q r + n r ■ qj = Vj, r G k (2-1) 

then there exists a reflected Brownian motion which can be defined as the unique solution 
to a submartingale problem. It was also shown in [36] that, under the skew-symmetry 
condition, reflected Brownian motion has an an invariant measure in product form: 

Theorem 2.2. Suppose that the skew-symmetry condition (2.1) holds, then RBM cor- 
responding to (N, Q, fi, b) has a unique invariant measure whose density with respect to 
Lebesgue measure is given by 

ir(x) = exp{-27(/i) ■ x} (2.3) 

where 7(/u) is defined as follows. By assumption, N has an invertible d x d submatrix 
N. Denote the corresponding submatrix of Q by Q, then 

1 ( Li ) = (l-N- 1 Qy 1 f ,. (2.4) 



3 



2 Reflected Processes in a Polyhedral Domain 



Remark 2.5. The existence of an invertible submatrix N of N was assumed. By the 
remarks after equation (1.7) in [16] (p. 463) the matrix (I — N Q), and hence 7(/x), 
is independent of the choice of N, provided the skew-symmetry condition (4.1) holds. 
Further [16, (4.7), (7.13)] we have \ j{fi)\ 2 = 7 • \i. 

Harrison- Williams [16] also consider reflected Brownian motion in a smooth 
domain and establish similar results in that setting. 

Remark 2.6. In two dimensions, the skew-symmetry condition corresponds to a con- 
stant angle of reflection across the entire boundary. To be more precise (cf. [16], Example 
8.1), fix an orientation on dG (which, in two dimensions, is always connected), and let To- 
be the unit tangent to the face Fa pointing in the positive direction with respect to this 
orientation. Define the angles 9j G (—§,§) by qj = tj tan#j. Then the skew-symmetry 
condition is equivalent to the requirement that 9j = 9 r for any r,j. 



2.2 RBM in an Orthant 

We now turn to reflected Brownian motion in the orthant S = (0,oo) a! , driven by 
a general-covariance Brownian motion. Our definition almost exactly mirrors that of 
Harrison-Reiman [15]. However, we have changed the sign of the reflection matrix 
Q to make it compatible with the Harrison- Williams setup. 

Let d G N and B be a ti-dimensional Brownian motion with drift [i and covariance 
matrix A = aa T ', started inside S. That is, there exists a A;- dimensional standard 
Brownian motion /3 and a k x d matrix a with unit rows such that B(t) = cr/3(t) — [it. 
In [15] additional assumptions on the matrix Q are required, namely that the dxd matrix 
Q has non-positive entries, spectral radius less than one and zeroes on the diagonals. 
Harrison-Reiman prove that under these assumptions there exists a unique pair of 
continuous M d -valued processes (Y, Z) with 




Figure 2.2: The skew-symmetry condition in two dimensions 



Z(t) = Y(t) + B(t)(I + Q) 



and such that 



(i) Z(t) e S for all t > 
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(ii) for each j G d the real-valued process Yj is continuous, non-decreasing and such 
that Yj(0) = 

(iii) Yj only increases at such times t where Zj(i) = 0. 

The process Z is called reflected Brownian motion in the orthant S with respect to the 
matrix Q, driven by B. 

If k = d and a is an invertible matrix then it is easy to see that the process a -1 (Z) is a 
reflected Brownian motion in the polyhedral domain a~ 1 (S) in the sense of Harrison- 
Williams. 



3 An Analytic Result 

Before introducting the GRBM process in the next section we state a purely analytic 
result. We continue to use the same notation and assumptions for the polyhedral domain 
G and the normal and reflection matrices N and Q. As mentioned in the introduction, 
the idea is to replace the singular reflection drift by a continuous drift that depends, via 
a function U on the position of the process relative to the domain G. This corresponds 
to defining a differential operator C by 

1 k 

£= =A-fi- V + ^2u'(n r -x-b r )v r ■ V (3.1) 

2 r=l 

and defining the GRBM process as diffusion with generator C. Let now v be a measure 
on R rf which is absolutely continuous with respect to Lebesgue measure and has density 
p and such that C*p = 0, where C* denotes the formal adjoint of C. That is, C* = 
^ A - ft • V - V ft where 

k 

fl(x) = U' (n r • x — b r ) v r — n 

r=l 

If the operator C is sufficiently well-behaved (which will depend on the function U) 
this would imply that the measure v is invariant with respect to the GRBM process. 
We begin by the following analytic result, which states that under the skew-symmetry 
condition (2.1) of HARRISON- WILLIAMS, the condition C*p = holds for a particular 
density function that can be considered the analogue of the function 7r from Section 2. 

We remark that the result does not require U to be anything more than twice 
continuously differentiable. In the next section we will discuss further conditions on U 
that ensure a probabilistic interpretation. The partial differential equation C* f = can 
be viewed as the analogue of the basic adjoint relations (BAR) of [16,36]. 

Theorem 3.2. Suppose that the nj,qj satisfy the skew-symmetry condition (2.1), let 
U : R — > R be any twice continuously differentiable function and define the function 
p:R d — > R by 



Pu 



(x) = exp 1 2 U{n r -x- b r ) - ■ x^j j dx (3.3) 



where j(fi) is as in (2.4)- Then C*p = 0. 
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The proof of this result is deferred to Section 6. A similar result holds for the analogue 
of RBM in an orthant as defined by Harrison-Reiman [15] (see Section 2.2). Here 
we allow the driving Brownian motion to have a more general covariance matrix, but 
require that rij = ej, where e\, . . . , ea denotes the standard orthonormal basis for H d . 
Therefore N = I, the d x d identity matrix. 

Theorem 3.4. Let A = (ajk)j t k be a symmetric positive definite matrix whose diago- 
nal entries are ajj = 1 for any j. Suppose that qi,---,qa satisfy the modified skew- 
symmetry condition 

Then we have C*p(x) = where C* denotes the formal adjoint of the operator 



(3.5) 



1 k 
C = - V • (AV) - y. ■ V + u ' ( x r) (Qr + e r ) • V 



(3.6) 



r=l 



and the density function p is defined by 

d 



p(x) = exp < 2 



Y J U(x J )-(2A-I-Q)- 1 f i 



> . 



(3.7) 



We note that both the original and the modified skew-symmetry condition do not depend 
on the function U. 



4 Generalised Reflected Brownian Motion 
4.1 The GRBM process 

We now turn to the probabilistic interpretation of the results from Section 3. Let B 
be a Brownian motion with drift —fi and covariance matrix A = (a,jk)j,k = O0 ~ ■> whose 
diagonal entries are cijj = a > for any j. The generator of B is given by \ V • (AV). 
We always tacitly assume that either the covariance matrix A = I (the Harrison- 
Williams setting), or d = k and N = I (i.e. the domain is an orthant, generalising 
the Harrison-Reiman RBM). We also assume that the generalised skew-symmetry 
condition 

rij ■ q r + n r ■ qj = — ^ Vr ^ j. (4-1) 

holds. Of course equation (4.1) reduces to (2.1) and (3.5) in the Harrison-Williams 
and Harrison-Reiman cases respectively. In order to interpret Theorems 3.2 and 3.4 
probabilistically we need the function U to satisfy a certain regularity condition. 

Definition 4.2. A twice continuously differentiable function U : M — > M is said to be 
regular with respect to the input data (N, Q,fi,A) if 

(1) the second-order differential operator C defined by 

1 k 
£ = - V ■ (AV) - ii ■ V + ^ U' (n r ■ x - b r ) v r ■ V (4.3) 

r=l 

is the infinitesimal generator of a diffusion process with continuous sample paths in 

R d 
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(2) if p: R d — > [0,oo) is smooth, integrable and such that C*p = 0, then p(x) dx is an 
invariant measure for this diffusion. Here C* denotes the formal adjoint of C. 

Definition 4.4. Suppose that the function U is regular with respect to the input 
data (N, Q, p,A). Then the Revalued diffusion with infinitesimal generator C is called 
generalised reflected Brownian motion ( GRBM) corresponding to the potential U and 
the data (N, Q, p, A). 

We now state a sufficient condition U to be regular with respect to (N,Q, p,, A) 
where one of TV or A are equal to the identity-covariance case. The proof is deferred to 
Section 6. 

Proposition 4.5. Suppose either A = I or k = d and N = I, that the nj,qj are as 
above, in particular satisfying the skew- symmetry condition (4-1) and that U satisfies 
the following conditions: 

(i) x — U(x) — > 00 as \x\ — > 00 

(ii) For each a£l there exists a > such that U{x + a) — U(a) < 9 a x for all x 6t. 

(Hi) There exist k, a > such that j r U'(x) — \ U"{x) < a {6-b r x — U(x)) + k for all 
x£R and r G k. Here, j r = Yl^=i ®-b s ( v r " n s) + p ■ n r . 

Then U is regular with respect to (N,Q, p,A). 

If N, Q satisfy the skew-symmetry condition then the choice U(x) = —e~ x satisfies the 
hypotheses of Proposition 4.5, and hence is regular with respect to (N,Q, p, A). We 
refer to the process with this choice of U as the exponentially reflected Brownian motion 
(ERBM) corresponding to (N,Q,p,A). 

Remark 4.6. Apart from the requirement that the rij contains a basis of M d one 
could allow k < d, that is fewer half-spaces than the dimension we are in. We can 
then decouple the 'superfluous' dimensions as follows. Denote by E\ = R h the span of 
{ni, . . . , rifc} and by E2 its orthogonal complement in M d . Let further Pj be orthogonal 
projection from M. d onto Ej. By orthogonal invariance of Brownian motion it follows 
that P\(X) and are independent. Further Pi(X) is just a standard (d — k)- 

dimensional Brownian motion, whereas P2(X) is generalised RBM in R k with the data 
(Pi(JV),Pi(Q),Pi(/i),/). 

With this in mind we will assume throughout that k > d. 

4.2 Invariant measures for GRBM in a general domain 

Having defined the GRBM process let us now establish the probabilistic interpretations 
of Theorems 3.2 and 3.4. We first consider the identity-covariance case A = I. The 
main result is that, under the skew-symmetry condition, generalised reflected Brownian 
motion has an invariant measure in a certain product form. It follows directly from 
Theorems 3.2 

Corollary 4.7. Suppose that U is regular with respect to (N, Q, p, I) and that the nj,qj 
satisfy the skew-symmetry condition (4-1)- Then the ~M. d -valued diffusion with generator 
C as defined in (3.1) has as invariant measure of the form vu{dx) = p(x)dx, where p 
is as in (3.7). 
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Remark 4.8. Under normal reflection (i.e. Q = 0) the GRBM process is a gradient 
diffusion, that is its generator is of the form \ A + VW ■ V. Note that Q = implies 
7(/u) = fi. So in this case our result also follows from the well-known fact (see [12] for 
example) that a gradient diffusion has e 2W ^ dx as invariant measure. 

Remark 4.9. If the function U satisfies the condition 

- l -U"{x) + {l + ^)U'{xf-^<- 1 

for all x £ and the corresponding density p is integrable then the stationary distri- 
bution obtained by suitably normalising (3.3) is actually unique. For example if k = d 
and U (x) = —e~ x , then there is a unique invariant measure if and only if fx ■ rij > for 
all j. 

Applied to the special case U(x) = —e~ x , Theorem 4.7 gives the invariant measure for 
exponentially reflected Brownian motion. 

Corollary 4.10. Suppose that the skew- symmetry condition (4-1) holds and the measure 



v(dx) = exp < 




> dx 



is finite. Then the exponentially reflecting Brownian motion corresponding to (N, Q, fi,I) 
has an invariant measure v. 

Remark 4.11. When d = k then v can be realised as the distribution of a M^-valued 
random variable \ such that 

( ni -x,---,n d -x) = ( _1 °s(y) '•••' _lo s(y)) 

where the jj are independent gamma random variables with parameters Oj and the 
vector 6 is given by 

= 2 (N T - Q T y X /i. 
That is, the law of nj ■ x is given by A^, where for a > 0, 

A a (da;) = \- exp {-2 [ax + e~ x ) } dx 

Remark 4.12. Let j3 > and set Up{x) = -^e~^ x for (3 > 0. The diffusion with 
generator C corresponding to U = Up should converge, as j3 — > oo to the Harrison- 
Williams reflected Brownian motion. Moreover (cf. [5], section 4.1) the log-gamma 
random variables converge to the exponential distribution, and we recover the main 
result of [36]. In this sense our results can be considered as a generalisation of those 
of [16,36]. 
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4.3 Invariant measure for GRBM in an orthant 



We now turn to exponential Brownian motion in an orthant, driven by a c?-Brownian 
motion with drift fi that is allowed to have a general, possibly singular covariance. 
It can be realised as B{t) = crf3(t) for a standard Brownian motion (5, possibly of a 
different dimension, and a matrix a, which is generally rectangular. The covariance 
matrix A = (ajk) = cra T of B has 1 on all its diagonal entries (or, equivalently, the 
rows of the rectangular matrix a have unit length) . This requirement can be weakened 
slightly, see Remark 4.16. 

In this case the generator of the generalised RBM is given by 



V • (AV) + 



d 

E 

r=l 



U' (x r ) (q r + e r ) - n 



V. 



The probabilistic interpretation of Theorem 3.4 is the following 

Corollary 4.13. Suppose that U : R — > R is regular with respect to (I,Q,n,A) and 
that for all j / r the modified skew-symmetry condition 



(£j rp | Q^fj 2 Ob ■ 



rj 



Vr / j 



(4.14) 



holds. Then the GRBM in an orthant corresponding to U has an invariant measure 
whose density with respect to Lebesgue measure is given by 



pu(x) = exp < 



d 

E 



U (xj) - (2A - I - Q)~ V ■ x 



(4.15) 



provided that f pu < oo. 

Remark 4.16. If the rows of A are not unit vectors but all have the same Euclidean 
norm ^Ja then we can scale the whole system so that our reasoning still applies. In this 
setting we get the stationary density in product form 



p(x) = exp < 2 



d 

E 



U 



a 



A-{I + Q)\ n-x 



(4.17) 



5 Examples 

5.1 One-dimensional ERBM and Dufresne's identity 

As a warm-up let us consider one-dimensional exponentially RBM. Here we can use a 
particular realisation of the process and Dufresne's identity. 

In this situation n = 1 and q = 0. All conditions, including skew-symmetry (4.1), 
are satisfied. Let further fi > 0. The generator of X in this simple case is given by 

4° = ^ («--") 5- 

By Ito's formula and stochastic integration by parts [17, 29] the process X given by 

X(t) = \og f e^'W-^C) ds 
Jo 
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is a diffusion with generator Cf . The invariant measure of X is that of n = log(£) 
where £ = and the process AW is defined by [7, 8] 

Recall Dufresne's identity [7, Corollary 4]. 
Proposition 5.1. Let fi> 0, then 

(2^) _1 ± 7, 
where 7 M /ias the Gamma distribution with parameter fi. 

Hence £ = ^- and so the invariant distribution of the process X is realised by 

v = -logC 1 = - log (Jy) ■ 

I.e. our results recover Dufresne's identity. Let us also note that, replacing the po- 
tential e~ x by ^e~^ x (as in Remark 4.12) it can be verified directly that both the 
one-dimensional reflecting Brownian motion and the exponential distribution appear in 
the scaling limit as (3 — > oo. 



5.2 TASEP-like particle systems 

Consider the following particle system on the line. There are n particles, with positions 
Xx(t), . . .,X n (t) at time t, labelled such that Xi(0) < X 2 (0) <■■•< X n (0). The left- 
most particle X\ evolves like a Brownian motion with drift u±, the second particle X 2 
like a Brownian motion with drift v<i reflected off X\ and so on. This particle system, 
first introduced by Glynn- Whitt [11], can be thought of as a Brownian analogue of 
the totally asymmetric simple exclusion process (TASEP). 

Replacing the 'hard' reflection by a 'soft' one, determined by the potential U we 
obtain a diffusion on W 1 which can be defined as the unique strong solution of the 
system of SDEs 

dXi(t) = dBi(t) 

dX k+1 (t) = dB k+1 (t) - U (X k+1 (t) - X k {t)) At (5.2) 

where B = {B\, . . . , B n ) is an n-dimensional Brownian motion with drift v £ W 1 . The 
choice of U{x) = —e~ x was first considered in [25]. By viewing the process as an 
analogue of a queueing system and applying a suitably generalised version of Burke's 
theorem a stationary distribution for the differences Yj = Xj + \ — Xj is obtained. In 
particular, if the drift of the driving Brownian motion is given by v = (a, 0, . . . , 0) then 
the stationary distribution of Y is given by ®j=i ^a- By ^ ne ^ aw °f l &r g e numbers, it 
follows that the speed of the particle system started in equilibrium is given by 

(a) = lim - X n (n) = Ed = - (5.3) 
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where T is the Gamma function and T' its derivative, the digamma function. This 
allows one to compute the free energy of a related semi-discrete directed polymer model. 
Suppose now that the drift of the driving Brownian motion B is zero and define 

Z n (t) = / exp{B 1 (s 1 ) + B 2 (s2)-B 2 (s 1 ) + ... + B n (t)-B n (s n - 1 )} dsi...ds n _i 

■Mn(t) 

where A n (t) = {(si, . . . , s n _i) G [0, oo) : s\ < . . . < s n _i < t}. It was shown in [25] 
that, almost surely, 

7 (a) := - (-*)* {a) = lim - log Z n (an) 

where ■* denotes convex conjugation, i.e. f*(a) = sup{ax — f(x): x G M.}. Note 
that (logZfc)fc=i satisfies (5.2) for the exponential case, with a particular entrance law. 
Hence, if we start the particle system X k with this entrance law (which corresponds in 
an appropriate sense to the 'step initial condition' X\(Q) = and Xj + i(0) — X/(0) = — oo 
for each j) we can compute the speed of the system: almost surely, 

lim — X n (na) = 7(a). (5.4) 

n— >oo n 

For a general potential U, the gaps between the particles, defined by Yj(t) = Yj(t) — 
Yj + i(t) evolve as a GRBM in the n — 1-dimensional orthant; hence, by Corollary 4.13, 
Y has an invariant measure given by 

^ n—l 

&/(dx) = — Y[ exp{-2 (U(x k ) + (y k - v x ) x k )} dx k 
u k=i 

where we assume that there exists Zjj G (0, 00) such that £u is a probability measure. 
So then it is still true that Y has a stationary distribution in product form, and under 
the special drift v = (a, 0, . . . , 0) the marginals are identically distributed. Hence the 
mean particle density of the system started in equilibrium is given, almost surely, by 

1 1 f 

\I/ U (a)= lim -X n (n) = — I x exp {-2U(x) + 2ax} dx. (5.5) 

n^oo n Zu J R 

provided that this is finite. Hence, under reasonable assumptions on U, it should be the 
case that, from an appropriate entrance law corresponding to the 'step initial condition' 
Xi(0) = and X k+ i(0) — X k (0) = —00 for each k, the speed of the system will be given 
by -HM*. 

5.3 Examples related to the Pitman transform 

The special case of U(x) = —e~ x in Section 5.2 is also related to the exponential 
analogue of the Pitman transform [1,2]. This transform is defined, for a G M ra and 
77 G C ((0, 00); M n ) by 

T aV (t) = rj(t) + (log e -« v -"W da) a, (5.6) 
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where a v = 2a/\a\ 2 . Suppose now that i] is a n-dimensional Brownian motion. For a 
sequence of vectors 71 , . . . , j q € W 1 and j € { 1 , . . . , q} define 

m = T 7 . o ■ ■ ■ o T 71 fa) = T 7j _ 1 (^-.x) (5.7) 

and then set 



= log f e 7/'fe-iW-<B-iW) ds . (5. 8 ) 
■/ 



We will see below (Proposition 5.10) that the process (yi, ■■ ■ ,y q ) fits in the framework 
of GRBM and hence has a product form invariant measure. 

Of particular interest is the following special choice of 71, . . . , j q . Let (Xj = ej — e.j+i, 
j = 1, . . . ,n — 1, where ei, . . . , e n denotes the standard orthonormal basis of W 1 . Let 
S(n) denote the symmetric group on n elements and denote the adjacent transpositions 
by Sj = (j j + 1) € S(n). Under the natural action of S(n) on R n , the Sj correspond 
to reflections in the orthogonal hyperplane to ctj. Let 



1 2 ... n 
n n — 1 ... 1 



be the longest element in S (n) , which has the following reduced decomposition of length 

n(n— 1) 

cr = si . . . s n _i si . . . s n _ 2 • • • sis 2 si- 
Let 71 , . . . , 7 g € R n be defined by 

71 . . . 7 ? = ai . . . a n -\ a\... a n -2 ■ ■ ■ ai«2 ol\ (5.9) 

and define the process y = (y±, . . . , y q ) as in (5.7) and (5.8). Then y has an invariant 
measure with density 



exp 



provided that /U • Q r > for all r. Here 0(/x) is a parameter which will be defined in 
Proposition 5.10, see (5.12). 

It follows from the definition of the yi that the stochastic evolution of the process 
(yi, ... , y n -i) is identical in law to the inter particle distances in the TASEP-like process 
X from Section 5.2 given the appropriate relation between the drifts \i and v. On the 
other hand, for any a £ S(n) and any reduced decomposition a = Sj 1 . . . Sj r the operator 
T a := T ajr o . . . o T aj only depends on a, not the choice of decomposition [1]. It was 
shown in [22] that the process ((T ao r]) (t)) t>0 is a diffusion which is closely related to 
the quantum Toda lattice, and can be thought of as a positive-temperature analogue 
of Dyson's Brownian motion. This allows one to compute the Laplace transform of 
the partition function of the semi-discrete directed polymer mentioned above. The 
process (yi, . . . ,y q ) can also be interpreted in terms of a natural diffusion on the set 
of totally positive lower triangular matrices [23]. The fact that this process has an 
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invariant measure in product form can be viewed as a positive temperature analogue of 
Proposition 5.9 in [2]. 

The appearance of product-form invariant measures in this context does not depend 
on the algebraic structure described above. In particular the above example has a 
natural generalisation which replaces S(n) with an arbitrary finite reflection group. We 
remark that a study of diffusions analogous to T ao (ry) for the generalised quantum Toda 
lattice (corresponding to arbitrary finite reflection groups) can be found in [4]. 

The following proposition describes the general setting. 

Proposition 5.10. Let j±, . . . , 7 g be any collection of vectors in W 71 such that \~/j\ 2 = 2 
and assume that U : R — > R is regular with respect to (I, T, /x, A) where Ajk = jj • 7^ 
and the matrix T is defined by 



r 



if j > k 
otherwise. 





Jj ■ Ik 

Suppose further that y = (y±, . . . , y q ) satisfies the system of SDEs 



'd( 7l - !?(*)) + [/' ( yi (t)) dt 

d (Tfc • m + feK (Tfc • 13) U' + U' (y k (t)) 



dt 



ifk = l 
ifk > 1 



(5.11) 



and i/iai 9 r • fi > for all r, where 6 r (fi) = —j3 r ■ fj,, the parameters (3 r are given by 



fa 

and the s 7 are defined by s 7 (x) - 
invariant measure with density 



71 



71 



-i(7fc) 



i/fc > 1. 



(5.12) 



-j^jf 7. T/ien i/te q-dimensional diffusion y has an 



> . 



(5.13) 



Proof. The process y is GRBM in a g-dimensional orthant, driven by the process n given 
by 



7i 



T 
v7g 



with /3 a standard Brownian motion. The covariance matrix of r/ is therefore given 
by A, while the reflection matrix is F. It is now straightforward to verify that the 
modified skew-symmetry condition holds. An application of Corollary 4.13 completes 
the proof. □ 

Even in one dimension, Proposition 5.10 yields non-trivial and interesting examples. 
For example, take U(x) = —e~ x and choose jj = (— 1) J+1 y/2. The corresponding 
system of SDEs for the yj is given by 

d yi (t) = V2dB{t) - e~ yi{t) dt (5.14) 

fc-i 

2j2(-l) j+k e- yj(t) -e- yi W dt. (5.15) 
j'=i 



dy k (t) = (-l) k+1 V2dB(t) + 
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A realisation of this sequence of one- dimensional processes is given by the following 
construction. Define the operator S on continuous real-valued paths on by 

(Srj) (t) = r)(t) - 2 log f ds. 

Jo 

Then we have 

(SB) (t) = B(t) - 2Jx(t) + 2Ji(0) 
where J\ is the exponential analogue of th future infimum in the Pitman theorem, i.e. 

/oo 
e~ B ^ du. 

If we define inductively sequences (Y n : n G N) and ( J n : n G N) of processes taking values 
in R by Y x (t) = SB(t) + Ji(t) and further 

/oo 
e S-B(u) du 

Y n+1 (t)=S n+1 B(t) + J n+1 (t) 

then Yi,Y 2 , . . . satisfies the system of SDEs (5.14). Thus, Ji(0), J2(0), J 3 (0), . . . is an 
i.i.d sequence of random variables and S acts on it by left shift. In particular, if the 
Brownian motion path (B(t), t > 0) could be recovered almost surely from the sequence 
«7i(0), J2(0), Js(0), . . ., then the system would be metrically isomorphic to a Bernoulli 
shift. We leave this as an open question. A discrete version of this open question was 
answered in the affirmative in [19]. 

5.4 Generalised Pitman transform for other functions 

The fact that invariant measures in product form appear irrespective of the choice of 
potential suggests generalising the exponential analogue of the Pitman transform to 
other functions than the exponential, as follows. 

Suppose that U: K — > R is such that, for smooth 77: (0, 00) — > R the differential 
equation 

d (V^r/)) (t) = dr](t) + U' (TW (rj) (t)) dt 

has a unique solution with T^ u \0+) = —00. This is true for U(x) = —e~ x , and for 
this choice of U we have = T^, where T a is the generalised Pitman transform 
as defined in (5.6). This is extended to higher dimensions by noting that the (expo- 
nential) generalised Pitman transfrom T a only acts on the linear span of a and leaves 
the orthogonal space a 1 - invariant. So we can define the generalised Pitman transform 
as follows. Write II a , Tl a ± for the orthogonal projections from M. d to the span and 
orthogonal complemement of a respectively and define by 

T^i 1 (t) = Ii a ^(t) + T^ (U aV (t))a. 
It could be interesting to study the algebraic properties of these operators, as in [1,2]. 
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5.5 RBM in a Weyl Chamber 

Reflected Brownian motion in the Weyl chamber 

n = jx G R d : x x > x 2 > ■ ■ ■ > , 

with normal reflection is a realisation of Brownian motion with rank- dependent drift, 
as studied for example by Pal-Pitman [26]. This is defined by taking d standard 
Brownian motions X\, . . . , Xj with increasing re-ordering Xn\, . . . , X/^ and drift \jl 
with — n £ Q such that at each time t the process Xq) has drift fij. It can also be 
viewed as a Doob transform of the Delta-Bose gas, see Prolhac-Spohn [27]. The 
corresponding GRBM has in fact been studied before, by Rost-Vares [30]. Since 
GRBM with normal reflection is a gradient diffusion it follows that this process has an 
invariant measure in product form as in Theorem 4.10 (see Remark 4.8). 



6 Proofs 



In this final section we prove Theorems 3.2 and 3.4 and Proposition 4.5. We first prove 
the purely analytic fact that the analogue of the basic adjoint relations hold in both 
the general domain (Theorem 3.2) and the orthant (Theorem 3.4) and then establish, 
by proving Proposition 4.5, a sufficient condition for U to be regular with respect to a 
given input data (N, Q, /j,, A). 

6.1 Proof of Theorem 3.2 

We wish to show that C*p = where C* is the formal adjoint of the generator C and 
p is as in (3.3) and (3.7) respectively. Recall from the comments after (3.1) that the 
formal adjoint of C is given by 

£* = -A-ft-V-V-ft 

2 



where the function ft : 



— > M is given by 

k 

Cl(x) = U' {n r • x — b r ) v r — \x. 



r=l 



Note further that p has the form p{x) = exp{W(x)} where 

k 



n r ■ x 



7 • x 



W(x) = 
Let us remark that 

n(s) = - VW(x) + (7 - n) + (nr ■ x - b 

3=1 



yj q r . 
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Because qj ■ rij = for all j we have V ■ fi = \ AW. Further Vp = pVW and therefore 
Ap = (AW + |VVF| 2 ) p. Hence it follows that 



1 k 

- VW + (7 - n) + U ' K ■ x - b r ) q r 



Vp - - Wp 



1 

2 



V- 2 |VVF| 2 p+(7-/i)-VVFp 



fc 1 
- ^ f7' (n r • x - 6 r ) <? r • Vp AWp. 



r=l 



Hence, 
C*p(x 



(/x - 7) ■ W(x) - Y, U K • » - W) q r ■ VW(x) 

p\. x ) r=1 

k 

= 2j2 u '( n j- x - b j) 0" - 7) ■ rij - 2 (/x - 7) ■ 7 



r=l 



— {/' (n r • x — 6 r ) C/' (n s • x — 6 S ) n r • q s — 2 [/' (n r • x — b r ) 7 • q s 

r,s 

k 

2j2u'(n r -x- b r ) [Nfi - (TV - Q) 7 ] r . 



r=l 



where we have used the skew-symmetry condition and [y] r denotes the r th entry of a 
vector y. The fact that the last line is equal to zero follows from the fact that Nfi = 
(N — Q) for any choice of invertible submatrix N of N (and corresponding submatrix Q 
of Q) and that each row of N must occur in at least one invertible submatrix, because 
of Lemma 6.1 below. 

In the proof above we have used the following simple lemma from linear algebra. 

Lemma 6.1. Let k > d and £ = {y±, . . . , y^} C M. d a set of vectors whose span is the 
whole ofM. d . If there exists r G k such that any collection of d elements of £ containing 
y r is linearly dependent then y r = 0. 

Proof. The assumption that £ spans M rf is equivalent to £ containing a basis. By re- 
ordering we may therefore assume that {yi, ■ ■ ■ ,Vd} is such a basis. If now the set 
{yi, . . . , yj, . . . , yd+i} (meaning the vector y r is left out) is linearly dependent for r £ d 
then yd+i is in the span of {yi, ■ • ■ ,y r , ■ ■ - Vd}- By assumption this is true for every 
r £ d, so 



Vd+i e Pi span {yi, . .. ,y r ,. . . ,y d } = {0} 
which completes the proof. 



r=l 



□ 
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6.2 Proof of Theorem 3.4 

Next we turn to the proof of Theorem 3.4, establishing the analogue of the basic adjoint 
relations for GRBM in an orthant. The generator of this process is now of the form 

1 



where ft is given by 



C= - V -AV + Q-V 



3=1 



The formal adjoint of C is given by C* = \ V • AV - Q • V - Vfi, and we need to prove 
that C*p = where p(x) = exp 



W(x) = 2^2u(x j )-S(jj) 
and 5(h) = (2A - I - Q)' 1 fi. Noting that 



1 d 
U(x) = - VW(x) + Y,U' ( Xj ) qj + 5(h) - h 



3=1 



and hence V ■ O = \ AW it follows that 



C*p =- (V ■ AV + VW ■ AVW) p 



P 



l VW + Y J U'{n r x)q\ VW 



3=1 



P 



+ 5(h)) ■ VWp - | AW. 



Dividing by p and then using VW(x) = 2 (^2j U' (rij ■ x) rij — 5(n)j we obtain 
= I V • (A - I) VW(x) + \ VW(x) -V(A-I) VW(x) 

pyx J —i Zi 

d 

-Y,U' (Xj) qj ■ VW(x) + (h~ 5(h)) ■ VW(x) 

3=1 

= 2 E [ u " te) ( A - %] + 2 E u ' (*;) ^ (*■■) ^ - /) ir 



- 2 U' (xj) [5(h) -(A -I) ej + e 3 -(A -I) 5(h)] + 25(h) -{A -I) 5(h) 

3=1 
d 

+ 2j2u'(x j )[q j -5(H) + (H j -S(H) j )} 

3=1 



2 U' (xj) U' (x r ) g jr + 25(h) ■ (S(ti) - h) 



3,r 
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The first term vanishes because the diagonal entries of A are all equal to 1. Using the 
fact that A is symmetric we have 



C*p(x) 
p{x) 



2j2u'( Xj )U'(x r ) \(A-I) jr -q jr 



3,r 



+ 2 £V (xj) [qj • + (jij - 5(rij) - 2e 3 -(A -I) 6(ji)] 

3=1 

+ 26(fj)-(A-I)S^) + 2S(n)-{S^)- fi ) 

The first line equals zero because of the skew-symmetry condition, whereas the other 
two lines vanish due to the definition of S(fi). 

This completes the proof of Corollary 4.13. □ 



6.3 Proof of Proposition 4.5 

Finally we prove the sufficient condition for U to be regular with respect to (N, Q,fi,A) 
stated in Proposition 4.5. 

By Section 4.3 in McKean [20] there exists a diffusion with that generator, up to 
a blow-up time r, and up to r the transition function of the diffusion has a smooth 
density. Assume that we have a C 2 -function V : R d — > R such that 

a) V is proper, i.e. there exist functions a, b: R — > R, strictly increasing to infinity 
and such that a(0) = 6(0) = and 

a(\x\) <V(x) <b(\x\) VxGM d 



b) there exist c > and e > such that CV < cV + e outside some compact set 

By Theorem 1 in Thygesen [35] (see also Gard [10], p. 131), this guarantees the 
existence of a unique diffusion X with generator C, which is a regular process, i.e. there 
is no blow-up in finite time, or r = oo almost surely. In particular the corresponding 
martingale problem is well-posed. Suppose now that p:R d — > [0, oo) is a smooth 
integrable function with C* p = 0. By Proposition 9.2 of Ethier-Kurtz [9] it is enough 
to show that 

J (Cf) (x)p(x)dx = 

for all / G Cc(R d ;R). But this follows from C*p = because of compact support and 
integration by parts. So the result holds as soon as we have a Lyapunov function V 
satisfying a) and b) above. We claim that the function V: R d — > [0, oo) defined by 

k 

V{x) = [0-br n r ■ x - (U (n r • x - b r ) - U (-b r ) )] 

r=l 

does the trick. Recall that we are only considering the cases where A = I or N = I. 
Suppose the former holds. We first verify b): using the skew-symmetry condition, the 
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fact that the rij are unit vectors and rij ■ qj = 0, 

CV(x) = - i Yj U " K • x - b r ) + J2 u ' {n r -x- b r ) 



r=l 
k 



r=l 



n r 



s=l 



U' (n r • x — b r ) rij 



r=l 



r=l 



< a ^2 [Q-b r {n r - x - b r ) — U(n r -x — b r )] - ^ 



n r 



r=l 



r=l 



< aV(x) + e 



for some constants a, e > 0, so that b) holds. To see that V is proper define a, b : [0, oo) — 
R by 

a(p) = inf V(x), b(p) = sup V(x). 

\A>p \x\< P 

Clearly the only point to prove is that a(p) — > oo as p — > oo, but this follows from 
the fact that the n r form a basis and y — U(y) — > oo as \y\ — > oo. So a) holds which 
completes the proof for A = I. 

Suppose now that k = d and N = I. Recall that the covariance matrix A of the 
driving Brownian motion has diagonal entries all equal to a > 0. By scaling we may 
assume without loss of generality that a = 1 . We now define a function V by 



So we have 



v(x) = J2[e oXr -(u(x r ) + u(o))]. 



r=l 



CV{x) = -\Y J U"{x j ) + Y J l3U'{x j ) 



-^2 U> ( X r) UrjU' (X r ) ~ ^ Pj0 
r,j j=l 
d r -, -i d 



i=i L 



1 



1jU'{x j )--U"{x j ) 



J2pj9 

3=1 



because A is a covariance matrix, hence non-negative definite. The same argument as 
for the case of A = I now shows that V satisfies a) and b). 
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